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ABSTRACT 

Wc have studied power density spectra (PDS) of 206 long Gamma- Ray Bursts (GRBs). 
We fitted the PDS with a simple power-law and extracted the exponent of the power- 
law (a) and the noise-crossing threshold frequency (/th)- We find that the distribution 
of the extracted a peaks around —1.4 and that of /th around 1 Hz. In addition, based 
on a sub-set of 58 bursts with known redshifts, we show that the redshift-corrected 
threshold frequency is positively correlated with the isotropic peak luminosity. The 
correlation coefficient is 0.57 ± 0.03. 
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1 INTRODUCTION 

Gamma-ray Bursts (GRBs) show very complicated time pro- 
files and, despite extensive investigations, are still not fully 
understood. The Fourier power density spectrum (PDS) of 
GRBs, on t he other hand, see m to show relatively sim- 
ple behavior. iGiblin et al.l (| 19981 ') found that a typical PDS 
shows a low-frequency power-law component and a high- 
frequency flat component (u s ually assoc iated with Pois- 
son noise). iBeloborodov et all (jl998l . l2000l) considered each 
GRB as a realization of some common stochastic process 
and showed that when averaged over many bursts the re- 
sulting PDS exhibits a power-law behavior with an exponent 
of -1.67, which the authors note is consistent with the -5/3 
Kolmogorov spectral index expected from processes involv- 
ing turbulent flow. In addition, they claim that there is a 
break in the averaged PDS at ~ 1 Hz. The authors were not 
in a position to correct their sample for the time dilation 
due to the c o smolo gical redshifts. 

lLazzatil (|2002t ) analyzed GRB power spectra by di- 
viding them into six lum inosity bins using the vari a bility- 
luminosity correlation (Fcnimorc & Ramirez- Ruia 120001 : 



Reichart et all l200ll: iGuidorzil 120051: iGuidorzi et all 12005.. 
20061 : iLi fe Paczvnskill2006l : IRizzuto et al.ll2007f l. The PDS 

was averaged in each bin after correcting for pseudo- 
redshift s obtained through the v ariability -lumino si ty re - 
lation (|Fenimore fc Ramirez-Ruiz! |2000| ). lLazzatil [|2002l ) 
showed that the dominant frequency (fd) of the PDS 
is strongly correlated with the variability parameter ob- 
tained by taking a modified variance of the de-trended light 
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curve (IFenimore fc Ramirez- Ruizll2000l ). Here f d is obtained 
by finding the maximum of the function / x PDS(/) (see 
lLazzatil (|2002l ) for more details). The author further states 
that the red-noise component of the averaged PDS for the six 
luminosity bins is well described by a broken power-law func- 
tion with a low- frequency slope of -2/3 and high-frequency 
slope of -2. In this case, the break frequency is a function of 
both luminosity and the variability parameter. 

iBorgonovo et all (|2007l ) did a similar analysis of power 
spectra but used measured redshift information to correct 
for time dilation effects before averaging. The burst sample 
was subdivided into two populations based on the calculated 
values of the autocorrelation function. After averaging, the 
PDS of one population shows a power-law index of ~ -2.0 
(consistent with the spectral index expected of Brownian 
motion) and the PDS of the other population is character- 
ized by a low-frequency exponentially decaying component 
and a high-frequency power-law component with an index of 
~ -1.6 (which again is consistent with the -5/3 Kolmogorov 
spectral index). 

Most of the previous work on power density spec- 
tra of GRBs has been based on observations with the 
Burst and Transient Source Experiment ( BATSE) on 
the C ompton Gamma R ay O bs ervat or y (|Giblin et al.l 
19981: iBeloborqdov et all Il998l . l200d : lLazzatil 12001 



Borgonovo et al.1 20071 ) where a relatively modest amount of 



redshift i nformation is availa ble. The launch of the Swift 
satellite (|Gehrels et al.l 12004 ) ushered in a new era of GRB 
research. Due to its rapidly disseminated, arcsecond GRB 
positions, Swift has enabled more subsequent redshift 
measurements of GRBs than ever before. The availability 
of redshift information enables the study of rest-frame 
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properties of bursts and provides an opportunity for further 
exploration of correlations involving burst parameters such 
as luminosity and variability. 

In this paper we present a study of Fourier power den- 
sity spectra of 206 Swift long bursts. Unlike previous work, 
we avoid averaging PDS of multiple bursts and examine 
them individually. We have developed a method to estimate 
the uncertainties in PDS for each burst. Then we extract 
PDS for all the GRBs in the sample and investigate the dis- 
tribution of the extracted parameters. The structure of the 
paper is the following: In section [2] we discuss our method- 
ology for extracting the PDS. In section [3] we present our 
results for a sample of 206 long bursts and investigate vari- 
ous correlations between the extracted parameters. In addi- 
tion, we propose a new frequency-luminosity relation based 
on a sample of 58 GRBs with spectroscopically measured 
redshifts. In section[4] we discuss observational biases of our 
results. Finally, in section [S] we summarize our conclusions. 
In this work, we have adopted the standard values for the 
cosmological parameters: Om = 0.27, Q,l ~ 0.73 and the 
Hubble constant Ho is 70 (kms _1 )/Mpc. Throughout this 
paper, the quoted uncertainties are at the 68% confidence 
level, unless noted otherwise. 




log Frequency (Hz) 

Figure 1. Power Density Spectrum of GRB 081203A. The 
low frequency power-law component is referred to as the 'red- 
noise' component and the flat high frequency region is called 
the 'white-noise' component. The inset shows the light curve of 
GRB081203A. 



2 DATA ANALYSIS 

2.1 Light Curve Extraction 

Swift BAT is a highly sen sitive, coded aperture instru- 
ment ijBarthelmv et al.H2005l 'l. BAT uses the shadow pattern 
resulting from the coded mask to facilitate few arc-minute 
localization of gamma-ray sources. In order to generate back- 
ground subtracted light curves, we used a process called 
mask weighting. The mask weighting assigns a ray-traced 
shadow value for each individual event, which then enables 
the user to calculate light curves or spectra. 

We used the batmaskwtevt and batbinevt tasks in 
the Swift BAT FTOOLS to generate mask weighted, 
background-subtracted light curves in the BAT energy range 
15 — 200 keV. The light curves that are generated have 
rates that are measured in counts per second per detector 
(counts/sec/det). In addition, the above tools also generate 
uncertainties associated with the rates that are calculated 
by propagation of errors from raw counts (subject to Poisso- 
nian noise). For the BAT instrument one can potentially go 
down to the minimum time binning of ~ 0.1 ms. However, 
in this work, we used 1-ms time binned light curves. 



2.2 Fourier Analysis 

We calculate the Fourier transform, a/, of each GRB 
light curve, R(t) (measured in counts/sec/detector), us- 
ing a standard Fast Fouri e r Transfor m (FFT]Q algo- 
rithm jjenkins fc Wattslll969l ; |Pressll2002l ). We used a time 
segment of the burst light-curve where the total fluence is 
accumulated (i.e. start and end times corresponding to burst 
T100 which is calculated by the battblocks task). The PDS 



1 We used the FFT routine in the IDL (Interactive Data Lan- 
guage) data analysis package. 
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of each burst is calculated using Pf = aj aZ. The power spec- 
tra are not normalized nor are they averaged. In addition, we 
have employed logarithmic binning for our power spectra. 

This pr ocess of treating PDS individually is different 
fro m that of iBeloborodov et all (| 19981 . 120001 ) ; iLazzatil (|2002r) 
and iBoreonovo et al.l ~ 20071 ). as they used some averaging 
process to obtain the slope of the red-noise component of 
the power spectra. The wide variety of light curves exhib- 
ited by GRBs is potentially indicative of different emission 
and scattering processes that eventually shape the observed 
light curves and therefore we have avoided averaging power 
spectra so as not to compromise this valuable information. 

The uncertainties of the individual PDS are calculated 
as follows: For each burst, we simulate 100 light curves 
based on the original light curve (R^n) an d its uncertainty 

{ nrcal error \ • „ 
(-"bin )> Le - 



rpsimulatcd 



rji'cal j. nreal error 
J^bin + S X -"-bin 



(i) 



Here £ is a random number generated from a gaussian distri- 
bution with the mean equal to zero and the standard devia- 
tion equal to one. For each simulated light curve we calculate 
a PDS. Then we re-bin each PDS logarithmically. The un- 
certainties in the original PDS (obtained from the original 
light curve) are derived by taking the standard deviation of 
the 100 simulated PDS. 

The power spectra for the GRBs in the sample are fit- 
ted with the function depicted in equation [2] (see figure[T]for 
a typical fit). This function consists of a power- law compo- 
nent (to fit the low-frequency "red noise" component) and 
a constant component (to fit the flat high-frequency "white 
noise" component). 



logP(/) 



"(log / - 

logP w 
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for / < / th 
for / > /th- 
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Here / tri is the threshold frequency where the red-noise com- 
ponent intersects the white-noise component of the power 
density spectrum, and P w is the white-noise power density. 
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Figure 2. The peak photon flux versus the T100 duration. The 
final sample with a significant red noise component is shown in 
red boxes and the sample with no or weak red-noise component is 
shown in blue diamonds. Note that conservative 10% uncertain- 
ties were assumed for T100 values. 



3 RESULTS 

Out of 451 GRBs which triggered Swift BAT from 2004 
December 19 to 2009 December 31, we selected a sample of 
226 long GRBs that show a significant red-noise component 
above the fiat white-noise region. In figure [5] we represent 
the two samples (the sample with clear red noise compo- 
nent is shown in red boxes and the sample with no or weak 
red-noise component is shown in blue inverted triangles) in 
a peak-photon-flux versus TlOO-duration plot. For the most 
part, the bursts that do not show a clear red-noise compo- 
nent are generally either weak and/or short in duration. 

For the selected sample of 226 long GRBs we fitted the 
corresponding PDS with a simple power-law behavior given 
in eq uation^ using the nonlinear least squares routine MP- 
FIT (|Markwardtll2009l ). A typical fit is shown in figured] 
Out of the 226 GRBs in the sample, 20 bursts could not 
be fitted by a simple power-law. These GRBs were excluded 
from further analysis. For the final sample of 206 bursts, 
the distributions of the extracted slopes (a) and threshold 
frequencies (/th) are shown in figure [3] and figure [4] respec- 
tively. The distribution of slopes (a's) has a Gaussian-like 
shape and peaks around ~-L4 with a of about 0.6. The 
distribution of threshold frequencies (/th's) peaks around 1 
Hz and also shows a broad distribution. The distribution 
of the redshift corrected f t h (i-e. /th(l + z)), as depicted in 
the bottom panel of figure [4] shows a large dispersion and 
non-gaussian shape. In figure [5] we show a plot of a and 
/th; we see a very weak positive correlation (0.24T0.02) but 
we note at this stage of the analysis that f t h has not been 
corrected for noise contamination nor has it been corrected 
for redshift. 

There are 76 GRBs in our sample with measured red- 
shifts (spectroscopic or otherwise). For this sub-sample it is 
interesting to see whether the extracted parameters are red- 
shift dependent. Figure [6] shows a (top panel), fth (middle 
panel), and the redshift-corrected fth (bottom panel) as a 
function of redshift. Very weak correlations are observed be- 
tween a and redshift and also between fth and redshift. How- 



Figure 3. Histogram of the extracted slopes (a). The distribution 
shows a peak around —1.4 ± 0.6. 
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Figure 4. A histogram of extracted threshold frequencies (fth) is 
shown in the top panel. The histogram peaks around ~ 1 Hz. The 
bottom panel shows a histogram of redshift-corrected f t h for the 
subset of bursts with redshift measurements. Both distributions 
show a large dispersion. 
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Figure 5. The extracted slope, a as a function of the threshold 
frequency, f t h- A very weak positive correlation is observed. 
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Figure 6. The slope a as a function of redshift (top panel), the 
threshold frequency as a function of redshift (middle panel), and 
redshift-corrected f t ^ as a function of redshift (bottom panel). 
The top and the middle panels show very weak correlations while 
the bottom panel shows no significant correlation. Note that con- 
servative 10% uncertainties were assumed for redshift values. 



ever, no significant correlation is observed between / t h(l + z) 

and re dshift. 

In lUkwatta et all (|2009l ) we proposed a correlation be- 
tween the isotropic peak luminosity and the red-shift cor- 
rected fth based on 27 GRBs. To investigate this further 
with a larger sample we have selected a sample of 58 
bursts with spectroscopically measured redshifts and good 
spectral information. For this sample, we have calculated 
isotro pic peak luminosity as described in lUkwatta et alj 
(2010). Based on the availability of spectral information, 
we have divided the sample into three sub-samples: "Gold" , 
"Silver", and "Bronze". The "Gold" sample wi th 15 bursts 
have all Band spectral parameters measured l|Band et al.l 
1993). In the "Silver" sample (15 bursts), the E p has been 
determined by fitting a cutoff power-law[j (CPL) to spec- 
tra. These 15 bursts do not have the high-energy spec- 
tral index, j3, measured, so we used the mea n value of the 
BATSE P distribution, whi ch is -2.36 ±0.31 (jKaneko et all 
20061 : ISakamoto et~ai1l2009l ) . The "Bronze" sample, with 28 
bursts, does not have a measured E p . We have estimated 
it using the power-law index (r) of a simple power-law 
(PL) fit as described in Sak amoto et al.l ((2009). For these 



dN/dE ~ E a exp (-(2 + a)E/E p ) 
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Figure 7. Isotropic peak luminosity as a function of the redshift- 
corrected threshold frequency, /th(l + - z ). The parameters are cor- 
related with a correlation coefficient of 0.77±0.02 and the best-fit 
power-law yields an exponent of 1.27 ± 0.12. 



28 bursts, the low-energy spectral index, a, and the high- 
energy spectral index, /3, were not known, so we used the 
mean values of the BATSE a and /3 distribut ion, which are 
-0.87 ± 0.33 and -2 . 36 ± 0.31 respectively (|Kaneko et all 
120061 : ISakamoto et al.ll2009T l . 

The isotropic luminosity as a function of the redshift- 
corrected threshold frequency is shown in figure[7] In the fig- 
ure, the "Gold", "Silver", and "Bronze" samples are shown 
in red, blue, and green filled circles respectively. A clear pos- 
itive correlation can be seen in the figure. The Pearson's cor- 
relation coefficient is 0.77 ±0.02, where the uncertainty was 
obtained through a Monte Carlo simulation. The probabil- 
ity that the above correlation occurs due to random chance 
is ~ 4.5 x 10 -9 . Our best-fit is shown as a red dashed line 
in figure [7] yielding the following relation between Li so and 
fth- 



(•3) 



logLiso = (51.79 ±0.12) + (1.27 ±0.12) log(/ th (z + 1)). 

To compensate for the large scatter in the plot, the uncer- 
ta inties of the fit parameters are multiplied by a factor of 
^/V/ndf = a/1766/56 ~ 6.0 where ndf is the number of 
degrees of freedom. The blue dotted lines indicate the esti- 
mated la confidence level, which is obtained from the cumu- 
lative fraction of the residual distribution taken from 16% 
to 84%. 

Our result for the slope in figu re [7] is consistent with 
the value of 1.4 ± 0.2 obtained by lUkwatta etafl (|2009T) 
usi ng 27 GRBs. T his is encouraging because the results 
of lUkwatta et all (|2009l ) were obtained using non-mask- 
weighted event-by-event data instead of the mask-weighted 
data that we use in the current work. We also note that 
with the increase of the sample size by about a factor of two 
the correlation coefficient has increased from 0.69 ± 0.03 to 
0.77 ± 0.02. The correlation between frequency and lumi- 
nosity is clearly intriguing but there remain observational 

biases which we address in a later section. 

It has been reported previously (|Beloborodov et al.l 

2000) that the PDS slope is correlated with the burst bright- 
ness. In order to check our sample for this effect, we display, 
in figure [8] the slope (a) against brightness indicators: the 
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Figure 8. The peak photon flux as a function of alpha (top panel) 

and the fluence as a function of alpha (bottom panel). Very weak Figure 10. Histogram of the peak photon flux of the sample of 

negative correlations are observed in both cases. 58 bursts used to generate figure [3 
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Figure 9. The peak photon flux as a function of threshold fre- 
quency (top panel) and the photon fluence as a function of thresh- 
old frequency (bottom panel). No significant correlation is ob- 
served between fluence and / t h but a significant correlation is 
observed between peak photon flux and / t h- 

peak photon flux and the fluence. Very weak negative cor- 
relations are observed in both cases. 

The other extracted parameter, the noise crossing 
threshold frequency (/th) of the PDS, is also expected to 
depend on the brightness of the GRB. Presumably, the 
'red-noise' component of the PDS comes primarily from the 
GRB but the flat 'white-noise' component can in principle 
arise from the Poisson noise (intensity fluctuations) asso- 
ciated with the GRB and the natural background in the 
field-of-view of the detector. For distant and/or intrinsically 
weak bursts, noise unrelated to the burst may dominate the 
observed white- noise component, thereby overwhelming the 
red-noise part of the signal. This, in turn, would make the 
extraction of the threshold frequency brightness-dependent. 
In figure [9] we plot the peak photon flux and the photon 
fluence as a function of the threshold frequency. The red 
dashed line in the top panel of figure [5] is the best fit, given 
by equation|4l and blue dotted lines indicate a la confidence 
interval. 



log PPF = (0.58 ± 0.03) + (0.82 ± 0.04) log / th (4) 

Indeed, a positive correlation can be seen between / t h and 
the peak photon flux. However, no significant correlation is 
observed between / t h and fluence. We discus this important 
matter further in the next section. 



4 DISCUSSION 

It is conceivable that the proposed frequency-luminosity cor- 
relation is a direct result of the observed correlation between 
/th and the peak photon flux of the burst (see the top panel 
of figure [9]l . If this is the case, then for a statistically signif- 
icant sample of bursts with similar apparent brightness, we 
should not see a correlation between / t h and Li so . In order 
to select a sample of GRBs with similar apparent brightness 
we plot in figure [10] the peak photon flux distribution for 
the sample of 58 bursts used to investigate the frequency- 
luminosity correlation in figure[7] We see from figure fTOl that 
about half of the sample (28 GRBs) have a very similar peak 
photon flux (0.0 < log (peak photon flux) < 0.5). For this 
subset of bursts we plotted their peak photon flux and Li so 
as a function of fth and the results are shown in figure [Til In 
the top panel of figure [TT] it is clear that there is a significant 
correlation between f t h and Li so with correlation coefficient 
of 0.60 ± 0.06. This implies that the correlation observed in 
figure [7] (/th(l + z) - Li so correlation) is not entirely due to 
the correlation seen in the top panel of figure [9] (/th— peak 
photon flux correlation) . We now correct the / t h of this lim- 
ited sample (with similar apparent brightness) for redshift 
to see its effect. Plotted in the bottom panel of figure [TTJ 
are the redshift corrected data. We note that the correlation 
strength increases to a value of 0.78 ± 0.04, in part due to 
the natural correlation between redshift and Li so . 

In addition, we can approach the issue from the other 
direction, i.e., we select a subset of bursts with similar lumi- 
nosity and ask the question whether the correlation between 
peak photon flux and /th comes from the proposed / t h(l + z) 
- Liso correlation. In order to perform this test we selected 
a subset of bursts which have the roughly the same Li so val- 
ues (51.5 < log Liso < 52.5) and plotted their peak photon 
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Figure 11. The isotropic peak luminosity as a function of thresh- 
old frequency (with and without redshift correction) for a sample 
of bursts with narrow apparent brightness range. The red-dashed 
and blue-dotted lines are the same fit curves shown in figure [7] 
The size of the star is proportional to the redshift of the burst. 



flux as a function of fth- In figure 1121 we show the peak 
photon flux as a function of fth (top panel) and the red- 
shift corrected fth (bottom panel). There is clearly a strong 
correlation between the two parameters in both panels. It 
is interesting, however, that after the redshift correction, 
the correlation strength drops significantly. Accordingly, it 
would appear that the correlation between the peak photon 
flux and fth (top panel of figure [9| is not entirely due to 
the /th(l + z) - Lisa correlation (figure [7]) . Since the spec- 
tral power is proportional to the square of the flux and the 
PDS follows a f~ a behavior (see figure [3]), we expect to see 
a correlation between peak photon flux and fth- Hence, this 
correlation is in most part observational. 

Now we turn to the question of the dependance of the 
extracted threshold frequency on the noise- level of the burst. 
The obvious question is how to determine the noise-level 
for each burst. One way of defining the noise- level is the 
following: 



Noise Level ; 



Std.Dev.(DetrendedLC) 
Peak Count Rate 



x 100%. 



(5) 



The detrending of the light curve (LC) can be done in 
a number of ways and we adopted the following method. 
We generated two light curves of the same burst with two 
bin sizes. In order to produce the coarser binned light curve, 
we chose a time bin size that resulted in at least 100 points 
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Figure 12. Peak photon flux as a function of fth (top panel) 
and redshift corrected fth (bottom panel) for a sample bursts 
with roughly constant luminosity. The red dashed line in the top 
panel is the best fit line obtained in figure [9] and blue dotted lines 
indicate lcr confidence interval. 



in the burst duration (T100). The other light curve may 
have bin sizes that vary from 1 ms up to the coarser bin 
size. Clearly, with the different binning, the two light curves 
will have a different number of points. In order to properly 
detrend, we need to have the same number of points in the 
two light curves. We accomplish this by using a simple linear 
interpolation of the coarser binned light curve. The interpo- 
lated light curve is then subtracted from the finer binned 
light curve to generate the detrended light curve. 

Using equation[5]we extract a noise level for each burst. 
However, the extracted noise-level depends on the bin size 
used in the detrending process. This aspect needs to be ei- 
ther removed or accounted for before the noise level of all 
the bursts can be treated on an equal footing. 

The level of the flat white-noise region of the PDS does 
not depend on the bin size, i.e., for a given burst the white 
noise level is constant irrespective of the bin size, and for 
that matter so too are the extracted parameters a and /th- 
in order to remove the bin size dependence in the extraction 
of the noise-level, we modify equation [5] as follows. 



Noise Level = 



1 Std.Dev. (Detrended LC) 



v 7 ^ 



Peak Count Rate 



x 100%. 



(6) 



Here N is the number of data points in the finely binned light 
curve. Our tests indicate that the results given by equation[6] 
do not depend on the time bin size of the light curve and 
provide a robust measure of the noise-level of a given burst. 

In order to further investigate the dependance of fth on 
the noise-level, we performed additional tests. We simulated 
different noise levels by adding increasing amounts of Gaus- 
sian noise to a burst light curve (in this case GRB 050315). 
Then we extracted fth values for each setting of the noise- 
level. Our results, the extracted frequency values versus the 
noise-level, are shown in figure [13] as a log-log plot. The 
threshold frequency does indeed depend on the noise-level. 
However, there is a linear relationship between the logarith- 
mic values of the two quantities. This relation is important 
to know because it can be used to correct the extracted fth 
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Figure 13. The extracted threshold frequency as a function of 
the noise-level in a log-log scale. The threshold frequency displays 
a power-law dependance on the noise-level of the burst with a 
index (A) of -0.84 ± 0.08 for GRB 050315. The inset shows the 
time profile of the burst. 



Figure 15. The threshold frequency versus the noise-corrected 
threshold frequency. The two parameters show a strong correla- 
tion with correlation coefficient of 0.86. The dashed line indicates 
the equality line of the two parameters. 
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Figure 14. Distribution of noise-levels (top panel) and noise-level 
slope, A, (bottom panel) in the sample. 



values to some nominal noise-level that is common to all 
bursts in the sample. 

By performing the same test on the other bursts in our 
sample, we established that the relation between f t h and the 
noise- level depends on the profile of the burst, i.e., the slope 
(A) of the log- log plot is different for each burst. Shown in 
the bottom panel of figure Q3] is the distribution the slopes, 
A, obtained for our sample of 58 bursts used in the fth-Li so 
relation. Correspondingly, the noise-level (NL) distribution 
is shown in the top panel of figure 1141 This distribution 
shows a clear peak around the log value of -0.2 (NL ~ 0.6) 
while the A distribution peaks around the value of -0.8. 

We are now in a position to treat all the bursts in our 
sample on an equal footing and test whether the fth-Li so 
correlation, observed in figure [7] survives. The aim is to 
extract threshold frequencies which are consistent with a 
noise-level that is common to all the bursts in our sample. In 
order to accomplish this, we choose an arbitrary noise level 
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Figure 16. The noise-corrected, redshift-corrected threshold fre- 
quency versus noise level. 



of NL =1.0 (see figure [13)) and use the following relation to 
extract a corrected fth for each burst: 

l°g /th[NL = l] = l°g /th[NL=burst] — Aburst log(NL[ NL = burst] ) • (7) 

Here, /th[NL=burst] is the extracted threshold frequency for 
a given burst, Aburst is the noise-level slope corresponding 
to the same burst and NLt NL=burst ] is the burst noise-level 
determined by equation [5] The correction procedure is re- 
peated for each burst in our sample. To gauge the size of the 
correction, we plot in figure [151 (in a log-log scale) the cor- 
rected / tn values versus the uncorrected fth- We note that 
there is a strong correlation between the two parameters. 
This is a reflection of the clustering of the NL and A seen 
in figure 1141 We also plotted the NL as a function of the 
noise-corrected, redshift-corrected / t h in figure [TBI There is 
no correlation between these two parameters, thus giving us 
confidence in the noise correction procedure. 

We show in figure [17] the noise-corrected threshold fre- 
quency - luminosity relation. As is evident, the relation sur- 



8 Ukwatta et al. 









- • Gold Sample 

'. • Silver Sample 

- • Bronze Sample 








_ * 


Jk-i i ~5 h 














: J 
- -) 


/ 

/ T 

1 

i i . , , . 


Correlotion : 0.57±0.03 - 

i " 


2 


-1 

log Corrected Threshold Freque 


1 ; 

ncy (1+z) 



Figure 17. The noise-corrected, redshift-corrected threshold fre- 
quency versus isotropic peak luminosity. The correlation coeffi- 
cient between the two parameters is 0.57 ± 0.03. The solid line 
shows the best-fit power-law with a index of 1.67 ± 0.01. The 
dashed line shows the best- fit from figure [7] 



Table 1. Correlation coefficients 



Coefficient Type 


Correlation Coefficient 


Null Probability 


Pearson's r 


0.57±0.03 


1.42 X 10" 5 


Spearman's r B 


0.58±0.04 


1.72 x 10" 6 


Kendall's r 


0.43±0.03 


2.03 x 10" 6 



vives the noise correction albeit with a somewhat smaller 
correlation coefficient of 0.57 ± 0.03. Various correlation co- 
efficients of the relation are shown in Table [1] where the 
uncertainties were obtained through a Monte Carlo simula- 
tion. The null probability that the correlation occurs due to 
random chance is also given for each coefficient type. 

The new best-fit is shown as a solid line in figure [17] 
yielding the following relation between Lj so and /th[NL=i] : 

logL iso = (52.2 ±0.1) + (1.67 ± 0.01) log(/ t hO+ 1)). (8) 

The uncertainties in the fitted parameters are expressed with 
the factor of ^f X 2 /ndf = ^/ 1255/56 w 5.0. 

5 CONCLUSION 

In this paper we have analyzed PDS of 206 GRBs. We 
fitted each PDS with a simple power-law and determined 
the red-noise exponent and the threshold frequency where 
white noise begins. For a subset of GRBs, we extracted 
a frequency-luminosity relationship. For this sample, we 
treated all bursts on an equal footing by determining a com- 
mon noise level, thereby minimizing the potential observa- 
tional biases. We summarize the main results of our analysis 
as follows: 

• The distribution of the extracted a (slope of the red- 
noise component) values peaks around -1.4 and that of /th 
around 1 Hz. 

• The dispersion in the distribution of a is large and so 



the Kolmogorov index of -5/3 is accommodated by our anal- 
ysis. 

• The distribution of the redshift-corrected threshold fre- 
quency shows a large dispersion and is non-gaussian in 
shape. 

• Evidence is presented for a possible frequency- 
luminosity relationship, i.e., the redshift-corrected f t h is cor- 
related with the isotropic luminosity. The correlation coef- 
ficient is 0.57 ± 0.03 and the best-fit power-law has an in- 
dex of 1.67 ± 0.01. We appreciate that in reality there may 
be complicated underlying interrelationships involving peak 
photon flux, /th, and redshift and therefore the evidence for 
the frequency-luminosity relation should be considered ten- 
tative. 

• The proposed frequency-luminosity correlation, if con- 
firmed, may serve to provide a measure of the intrinsic vari- 
ability observed in GRBs. 
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